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Abstract 

This paper reviews recent results on hybrid inverse problems, which are also 
called coupled-physics inverse problems of multi-wave inverse problems. Inverse 
problems tend to be most useful in, e.g., medical and geophysical imaging, when 
they combine high contrast with high resolution. In some settings, a single modal- 
ity displays either high contrast or high resolution but not both. In favorable 
situations, physical effects couple one modality with high contrast with another 
modality with high resolution. The mathematical analysis of such couplings forms 
the class of hybrid inverse problems. 

Hybrid inverse problems typically involve two steps. In a first step, a well- 
posed problem involving the high-resolution low-contrast modality is solved from 
knowledge of boundary measurements. In a second step, a quantitative reconstruc- 
tion of the parameters of interest is performed from knowledge of the point-wise, 
internal, functionals of the parameters reconstructed during the first step. This 
paper reviews mathematical techniques that have been developed in recent years 
to address the second step. 

Mathematically, many hybrid inverse problems find interpretations in terms 
of linear and nonlinear (systems of) equations. In the analysis of such equations, 
one often needs to verify that qualitative properties of solutions to elliptic linear 
equations are satisfied, for instance the absence of any critical points. This paper 
reviews several methods to prove that such qualitative properties hold, including 
the method based on the construction of complex geometric optics solutions. 
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1 Introduction 

The success of most medical imaging modalities rests on their high, typically sub- 
millimeter, resolution. Computerized Tomography (CT), Magnetic Resonance Imaging 
(MRI), or Ultrasound Imaging (UI) are typical examples of such modalities. In some 
situations, these modalities fail to exhibit a sufficient contrast between different types 
of tissues, whereas other modalities, for example based on the optical, elastic, or electri- 
cal properties of these tissues, do display such high contrast. Unfortunately, the latter 
modalities, such as e.g.. Optical Tomography (OT), Electrical Impedance Tomography 
(EIT) or Elastographic Imaging (EI), involve a highly smoothing measurement operator 
and are thus typically low-resolution as stand-alone modalities. 

Hybrid inverse problems concern the combination of a high contrast modality with 
a high resolution modality. By combination, we mean the existence of a physical mech- 
anism that couples these two modalities. Several examples of physical couplings are re- 
viewed in section [2J A different strategy, consisting of fusing data acquired independently 
for two or more imaging modalities, is referred to as multi-modality imaging and is not 
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considered in this paper. Examples of possible physical couplings include: optics or elec- 
tromagnetism with ultrasound in Photo- Acoustic Tomography (PAT), Thermo-Acoustic 
Tomography (TAT) and in Ultrasound Modulated Optical Tomography (UMOT), also 
called Acousto-Optic Tomography (AOT); electrical currents with ultrasound in Ul- 
trasound Modulated Electrical Impedance Tomography (UMEIT), also called Electro- 
Acoustic Tomography (EAT); electrical currents with magnetic resonance in Magnetic 
Resonance EIT (MREIT) or Current Density Impedance Imaging (CDII); and elasticity 
with ultrasound in Transient Elastography (TE). Some hybrid modalities have been ex- 
plored experimentally whereas other hybrid modalities have not been tested yet. Some 
have received quite a bit of mathematical attention whereas other ones are less well un- 
derstood. While more references will be given throughout the review, we refer the reader 
at this point to the recent books [5l |53l |62] and their references for general information 
about practical and theoretical aspects of medical imaging. 

Reconstructions in hybrid inverse problems typically involve two steps. In a first 
step, an inverse problem involving the high-resolution-low-contrast modality needs to 
be solved. In PAT and TAT for instance, this corresponds to reconstructing the initial 
condition of a wave equation from available boundary measurements. In UMEIT and 
UMOT, this corresponds in an idealized setting to inverting a Fourier transform that 
is reminiscent of the reconstructions performed in MRI. In Transient Elastography, this 
essentially corresponds to solving an inverse scattering problem in a time-dependent 
wave equation. In this review, we assume that this first step has been performed. 

Our interest is in the second step of the procedure, which consists of reconstructing 
the coefficients that display high contrasts from the mappings obtained during the first 
step. These mappings involve internal functionals of the coefficients of interest. Typi- 
cally, if 7 is a coefficient of interest and u is the solution to a partial differential equation 
involving 7, then the internal "measurements" obtained in the first step take the form 
H[x) = '-f{x)u^{x) for j = 1,2 or H{x) = 7(x)|Vup(x) again for j = 1,2. 

Several questions can then be raised: are the coefficients, e.g. 7, uniquely charac- 
terized by the internal measurements H{x)l How stable are the reconstructions? If 
specific boundary conditions are prescribed at the boundary of the domain of interest, 
how do the answers to the above questions depend on such boundary conditions? The 
answers to these questions depend on the physical model of interest. However, there are 
important common features that we would like to present in this review. 

One such feature relates to the stability of the reconstructions. Loosely speaking, 
an inverse problem is well-posed, or at least not severely ill-posed, when singularities in 
the coefficients of interest propagate into singularities in the available data. The map 
reconstructed during step 1 provides local, point- wise, information about the coefficients. 
Singularities of the coefficient do not need to propagate to the domain's boundary and 
we thus expect resolution of hybrid modalities to be significantly improved compared 
to the stand-alone high-contrast-low-resolution modalities. This will be verified on the 
examples reviewed here. 

Another feature is the relationship between hybrid inverse problems and nonlinear 
partial differential equations. Typically, both the coefficient 7 and the solution u are 
unknown. However, for measurements of the form H[x) = '-f{x)u^{x), then 7 in the 
equation for u can be eliminated using the expression for H{x). This results in a 
nonlinear equation for u{x). The resulting nonlinear equations often do not display any 
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of the standard features that are amenable to proofs of uniqueness, such as admitting a 
variational formulation with a strictly convex functional. The main objective is to obtain 
uniqueness and stability results for such equations, often in the presence of redundant 
(overdetermined) information. 

A third feature shared by many hybrid inverse problems is that their solution strate- 
gies often require that the forward solution u satisfy certain qualitative properties, such 
as for instance the absence of any critical point (points where Vu = 0). The deriva- 
tion of qualitative properties such as lower bounds for the modulus of a gradient is a 
difficult problem. In two dimensions of space, the fact that critical points of elliptic 
solutions are necessarily isolated is of great help. In higher dimension, such results no 
longer hold in general. A framework to obtain the requested qualitative behavior of the 
elliptic solutions is based on the so-called complex geometric optics (CGO) solutions. 
Such solutions, when they can be constructed, essentially allow us to treat the unknown 
coefficients as perturbations of known operators, typically the Laplace operator. Us- 
ing these solutions, we can construct an open set of boundary conditions for which the 
requested property is guaranteed. This procedure provides a restricted class of bound- 
ary conditions for which the solutions to the hybrid inverse problems are shown to be 
uniquely and stably determined by the internal measurements. From a practical point 
of view, these mathematical results confirm the physical intuition that the coupling of 
high contrast and high resolution modalities indeed provides reconstructions that are 
robust with respect to errors in the measurements. 

The rest of this paper is structured as follows. Section |2] is devoted to the modeling 
of the hybrid inverse problems and the derivation of the internal measurements for the 
applications considered in this paper, namely: PAT, TAT, UMEIT, UMOT, TE, CDII. 
The following two sections present recent results of uniqueness and stability obtained for 
such hybrid inverse problems: Section 12] focuses on internal functionals of the solution 
u of the forward problem, whereas section S] is concerned with internal functionals of 
the gradient of the solution Vn. As we mentioned above, these uniqueness and stability 
results hinge on the forward solutions u to verify some qualitative properties. Section 
|5] summarizes some of these properties in the two-dimensional case and presents the 
derivation of such properties in higher spatial dimensions by means of complex geometric 
optics (CGO) solutions. Some concluding remarks are proposed in section El 

2 Physical modeling 

High resolution imaging modalities include Ultrasound Imaging and Magnetic Reso- 
nance Imaging. High contrast modalities include Optical Tomography, Electrical Impe- 
dance Tomography, and Elastography. This sections briefly presents four couplings 
between high-contrast and high-resolution modalities: two different methods to couple 
ultrasound and optics or (low frequeny) electromagnetism in PAT/TAT via the photo- 
acoustic effect and in UMOT/UMEIT via ultrasound modulation; the coupling between 
Ultrasound and Elastography in Transient Elastography; and the coupling between Elec- 
trical Impedance Tomography and Magnetic Resonance Imaging in CDII/MREIT. 
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2.1 The Photo- acoustic effect 



The photoacoustic effect may be described as follows. A pulse of radiation is sent into a 
domain of interest. A fraction of the propagating radiation is absorbed by the medium. 
This generates a thermal expansion, which is the source of ultrasonic waves. Ultrasound 
then propagates to the boundary of the domain where ultrasonic transducers measure 
the pressure field. The physical coupling between the absorbed radiation and the emitted 
sound is called the photoacoustic effect. This is the premise for the medical imaging 
technique Photoacoustic Tomography (PAT). 

Two types of radiation are typically considered. In Optoacoustic Tomography (OAT), 
near-infra-red photons, with wavelengths typically between QOOnm and 900nm are used. 
The reason for this frequency window is that they are not significantly absorbed by wa- 
ter molecules and thus can propagate relatively deep into tissues. OAT is often simply 
referred to as PAT and we will follow this convention here. In Thermoacoustic Tomog- 
raphy (TAT), low frequency microwaves, with wavelengths on the order of Im, are sent 
into the medium. The rationale for using such frequencies is that they are less absorbed 
than optical frequencies and thus propagate into deeper tissues. 

In both PAT and TAT, the first step of an inversion procedure is the reconstruction of 
the map of absorbed radiation from the ultrasonic measurements. In both applications, 
the inversion may be recast as the reconstruction of an initial condition of a wave 
equation from knowledge of ultrasound measurements. Assuming a domain of infinite 
extension with non-perturbative measurements to simplify the presentation, ultrasound 
propagation is modeled by the following wave equation 



Here Cs is the sound speed assumed to be known, n is spatial dimension, and H{x) is 
the ultrasonic signal generated at time t = 0. Measurements are then of the form p{t, x) 
for t > and x G dX at the boundary of a domain X where H{x) is supported. 

Note that the effect of propagating radiation is modeled as an initial condition at 
t = 0. The reason for this stems from the large difference between light speed (roughly 
2.310^m/s in water) and sound speed (roughly 1.510^m/s in water). When a short 
pulse of radiation is emitted into the medium, we may assume that it propagates into 
the medium at a time scale that is very short compared to that of ultrasound. This 
is a very valid approximation in PAT but is a limiting factor in the (still significantly 
sub-millimeter) spatial resolution we expect to obtain in TAT; see, e.g., [T2| [T6]. 

For additional references to the photo-acoustic effect, we refer the reader to the 
works [221 1231 [281 [SH [65] and their references. The first step in thermo- and photo- 
acoustics is the reconstruction of the absorbed radiation map H{x) from boundary 
acoustic wave measurements. There is a vast literature on this inverse source problem 
in the mathematical and physical literatures. We refer the reader to e.g. [H [271 [SH [33 
H2l [501 156] . Serious difficulties may need to be addressed in this first step, such as e.g. 
limited data, spatially varying acoustic sound speed [3 [351 [56], and the effects of acoustic 
wave attenuation [39]. In this paper, we assume that the absorbed radiation map H{x) 
has been reconstructed. This provides now internal information about the properties 
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of the domain of interest. What we can extract from such information depends on the 
model of radiation propagation. The resulting inverse problems are called Quantitative 
PAT (QPAT) and Quantitative TAT (QTAT) for the different modalities of radiation 
propagation, respectively. 

In the PAT setting with near-infra-red photons, arguably the most accurate model 
for radiation propagation is the radiative transfer equation. We shall not describe this 
model here and refer the reader to [12] for QPAT in this setting and to [10] for more 
general inverse problems for the radiative transfer equation. The models we consider 
for radiation propagation are as follows. 



2.1.1 QPAT modeling 

In the diffusive regime, photon (radiation) propagation is modeled by the following 
second-order elliptic equation 

—V ■ 7(x) Vn + cr(x)n = in X 

u = f on dX. ^ ^ 

To simplify, we assume that Dirichlet conditions are prescribed at the boundary of the 
domain dX. Throughout the paper, we assume that X is a bounded open domain in 
with smooth boundary dX. The optical coefficients {'y{x),a{x)) are 7(x) the diffusion 
coefficient and cr{x) the absorption coefficient, which are assumed to be bounded from 
above and below by positive constants. 

The information about the coefficients in QPAT takes the following form: 

H{x) = T{x)a{x)u{x) a.e. x G X. (3) 

The coefficient T{x) is the Griineisen coefficient. It models the strength of the photo- 
acoustic effect, which converts absorption of radiation into emission of ultrasound. The 
objective of QPAT is to reconstruct (7, a, T) from knowledge of H{x) in Q obtained for 
a given number of illuminations / in ([2]) . This is an example of an internal measurement 
that is linear in the solution u{x) and the absorption coefficient a. For references on 
QPAT, see, e.g., [H [TH |22l [231 Ell [66] and their references. 



2.1.2 QTAT modeling 

Low frequency Radiation in QTAT is modeled by the following system of Maxwell's 
equations: 

-VxVx E + k^E + tka{x)E = 0, X 

vxE = f dX. ^ ^ 

Here, E is the (time-harmonic) electromagnetic field with fixed wavenumber k = - 
where u is frequency and c light speed. We assume that radiation is controlled by the 
boundary condition /(x) on dX. The unknown coefficient is the conductivity (absorp- 
tion) coefficient cr(x). Setting F = 1 to simplify, the map of absorbed electromagnetic 
radiation is then of the form 

H{x) = a{x)\E\'^{x). (5) 
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The above system of equations may be simplified by modeling radiation by a scalar 
quantity u{x). In this setting, radiation is modeled by the following Helmholtz equation: 

Au + k'^u + ika{x)u = 0, X 

u = f dX, ^ ^ 

for a given boundary condition f{x). The internal data are then of the form 

H{x) =a{x)\uWx). (7) 

For such models, QTAT then consists of reconstructing a{x) from knowledge of H{x). 
Note that H{x) is now a quadratic quantity in the solutions E{x) or u{x). There are 
relatively few results on QTAT; see 



2.2 The ultrasound modulation effect 

We consider the following elliptic equation 

- V ■ 7(x)Vm + (t{x)u = in X, u = f on dX. (8) 

The objective of ultrasound modulation is to send an acoustic signal through the domain 
X that modifies the coefficients 7 and a. We assume here that the sound speed is 
constant and that we are able to generate an acoustic signal that takes the form of the 
plane wave pcos{k ■ x + (f) where p is amplitude, k wave-number and (p an additional 
phase. We assume that the acoustic signal modifies the properties of the diffusion 
equation and that the effect is small. The coefficients in ([8]) are thus modified as 

7.(x) = 7(x)(l + C^c) + 0(£2), a,ix) = a{x)il + r,ec) + ©(e^), (9) 

where we have defined c = c(x) = cos{k ■ x + ip) and where e = pF is the product 
of the acoustic amplitude p G M and a measure F > of the coupling between the 
acoustic signal and the modulations of the constitutive parameters in ([H]). We assume 
that e <C 1. The terms in the expansion are characterized by C, and 77 and depend on 
the specific application. 

Let u and v be solutions of ([8]) with fixed boundary conditions / and h, respectively. 
When the acoustic field is turned on, the coefficients are modified as described in ([9]) 
and we denote by and the corresponding solution. Note that is the solution 
obtained by changing the sign of p or equivalently by replacing by + vr. 

By standard regular perturbation arguments, we find that Ue = uq + eui + O(e^). 
Multiplying the equation for Ug by V-g and the equation for f_£ by Wg, subtracting the 
results, and using standard integrations by parts, we obtain that 

f f Ov-^ du 

I - 7-e)VM£ ■ VV-e + (o-e - (T-e)UeV-edx = / 7_e— -^M^ - '^^^V^edcr. (10) 

Jx Jdx du dv 

We assume that 'jed^Ue and 'jed^Vs are measured on dX, at least on the support ofv^ = h 
and Ue = f, respectively, for several values e of interest. The above equation also holds 
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if the Dirichlet boundary conditions are replaced by Neumann boundary conditions. Let 
us define 

Js :=]: [ - le^v-sda = eJ, + e'J, + 0{e'). (11) 

2 Jqx du du 

We assume that the real valued functions Jm = Jm{k,f) are known from the physical 
measurement of the Cauchy data of the form {u^.'^edyUe) and {v^,'~^^dyV^) on dX. 
Equating like powers of e, we find that at the leading order 

\C,'y{x)VuQ ■ Wvq{x) + ria{x)uQVQ{x)\ cos{k ■ x + ip)dx = Ji{k, (p). (12) 

Acquiring this for all k G M" and <y9 = 0, |, this yields after inverse Fourier transform: 

H[uo,vo]{x) = C-f{x)Vuo ■ Vvo{x) + r]a{x)uoVo{x). (13) 

In the setting of ultrasound modulated optical tomography (UMOT), the coefficients 
7e and in take the form [T7] 7e(x) = -^^{x) and o"e(x) = ^j|^(x), where is 
the absorption coefficient, % is the diffusion coefficient, is the light speed, and n 
is spatial dimension. When the pressure field is turned on, the amount of scatterers 
and absorbers is modified by compression and dilation. Since the diffusion coefficient is 
inversely proportional to the scattering coefficient, we find that 

11 

ae{x) = a{x){l+ec{x)), — - = —-{l+ec{x)). 

The pressure field changes the index of refraction of light as follows Cs{x) = c(x)(l + 
^£c{x)), where ^ is a constant (roughly equal to | for water). This shows that 

C = -(l + (n-l)7/;), 7] = l-{n-l)ij. (14) 

In the application of ultrasound modulated electrical impedance tomography (UMEIT), 
7(x) is a conductivity coefficient and a = 0. We then have 7e(x) = 7(x)(l + ec{x)) with 
thus C = 1 cind 7] = 0. The objective of UMOT and UMEIT is to reconstruct (part of) 
the coefficients (7(x),(t(x)) in the elliptic equation 

- V ■ 7(x)Vm + a{x)u = in X, u = f on dX. (15) 

from measurements of the form 

H[uo,vo]{x) = C-f{x)\/uo ■ Vvo{x) + r]a{x)uoVo{x), (16) 

for one or several values of the illumination f{x) on dX. 

In a simplified version of UMOT (also called acousto-optic tomography; AOT), C = 
and the measurements are quadratic (or bilinear) in the solutions to the elliptic equation. 
More challenging mathematically is the case C = 1 rj = where the measurements 
are quadratic (or bilinear) in the gradients of the solution. No theoretical results exist 
to date in the setting where both ( and t] are non-vanishing. 

The effect of ultrasound modulation is difficult to observe experimentally as the 
coupling coefficient F above is rather small. For references on ultrasound modulation in 
different contexts, we refer the reader to, e.g., [6l [9l [T71 [211 ISHl SH |67]. These references 
concern the so-called incoherent regime of wave propagation, while the coherent regime, 
whose mathematical structure is different, is addressed in the physical literature in, e.g., 

ilEZlEI]. 
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2.3 Transient Elastography 



Transient elastography images the (slow) propagation of shear waves using ultrasound. 
For more details, see, e.g., [H] and its extended list of references. As shear waves 
propagate, the resulting displacements can be imaged by ultra-fast ultrasound. Consider 
a scalar approximation of the equations of elasticity 

V ■-f{x)Vu{x,t) = p{x)dttu{x,t), teR,xeX 

u{x, t) = f{x, t), teR, X edX, ^ ' 

where u{x,t) is the (say, downward) displacement, 7(0;) is one of the Lame parameters 
and p{x) is density. Using ultra-fast ultrasound measurements, the displacement u{x,t) 
can be imaged. This results in a very simplified model of transient elastography where 
we aim to reconstruct (7, p) from knowledge of u{x, t); see [H] for more complex models. 
We may slightly generalize the model as follows. Upon taking Fourier transforms in the 
time domain and accounting for possible dispersive effects of the tissues, we obtain 

V ■ 7(x; u)'Vu{x] u) + uj'^p{x] uj)u{x] u) = 0, w e M, z e X 
u{x] u) = f{x] u), u eR, X e dX. 

The inverse transient elastography problem with dispersion effect would then be the 
reconstruction of {'j{x;u),p{x;u)) from knowledge of u{x;u) corresponding to one or 
several boundary conditions f{x; u) applied at the boundary dX. This hybrid inverse 
problem again involves measurements that are linear in the solution u. 



2.4 Current Density Imaging 

Magnetic Impedance Electrical Impedance Tomography (MREIT) and Current Density 
Impedance Imaging (CDII) are two modalities aiming to reconstruct the conductivity in 
an equation using magnetic resonance imaging (MRI). The electrical potential u solves 
the following elliptic equation 

- V ■ 7(x)Vm = in X, u = f on dX, (19) 

with 7(3;) the unknown conductivity and / a prescribed voltage at the domain's bound- 
ary. The electrical current density J = — 7Vn satisfies the system of Maxwell's equations 

V ■ J = 0, J = —V X B, xeX. (20) 

Here /iq is a constant, known, magnetic permeability. 

Ideally, the whole field B can be reconstructed from MRI measurements. This pro- 
vides access to the current density J{x) in the whole domain X. CDI then corresponds 
to reconstructing 7 from knowledge of J. In practice, acquiring B requires rotation of 
the domain of interest (or of the MRI apparatus) which is not straightforward. MREIT 
thus assumes knowledge of the third component B^ of the magnetic field for several pos- 
sible boundary conditions. This provides information about 7(0;). We do not consider 
the MREIT inverse problem further and refer the reader to the recent review [51] and 
its references for additional information. 
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Several works have considered the problem of the reconstruction of 7 in f ll9p from 
knowledge of the scalar information \J\ rather than the full current J. This inverse 
problem, referred to as the 1— Laplacian, will be addressed below and compared to the 
0— Laplacian that appears in UMEIT and UMOT. For references on MREIT and CDII, 
we refer the reader to SZl SHI SH] and their references. 

3 Reconstructions from functionals of u 

In this section, we consider internal measurements H{x) of the form H{x) = t{x)u{x) 
for t{x) a function that depends linearly on unknown coefficients such as the diffusion 
coefficient 7 or the absorption coefficient a in section 13.11 and internal measurements 
H{x) of the form H{x) = t{x)\u{x)\'^ in section [3l2| where r again depends linearly 
on unknown coefficients. Measurements of the first form find applications in Quan- 
titative Photo-acoustic Tomography (QPAT) and Transient Elastography (TE) while 
measurements of the second form find applications in Quantitative Thermo-Acoustic 
Tomography (QTAT) and simplified models of Acousto-Optics Tomography (AOT). 

3.1 Reconstructions from linear functionals in u 

Recall the elliptic model for photon propagation in tissues: 

- V ■ 7(x)Vn + a{x)u = in X, u = f on dX. (21) 

The information about the coefficients in QPAT takes the following form: 

H{x) = r{x)(T{x)u{x) a.e. x e X. (22) 

The coefficient r(x) is the Griineisen coefficient. In many works in QPAT, it is assumed 
to be constant. We assume here that it is Lipschitz continuous and bounded above and 
below by positive constants. 

3.1.1 Non-unique reconstruction of three coefficients. 

Let /i and /2 be two Dirichlet conditions on dX and ui and U2 be the corresponding 
solutions to fl?I]) . We make the following assumptions: 

(i) The coefficients (7, a, F) are of class W^'°°{X) and bounded above and below by 
positive constants. The coefficients (7, a, F) are known on dX. 

(ii) The illuminations /i and /2 are positive functions on dX and are the traces on 
dX of functions of class C^{X). 

(iii) the vector field 

(3 := - H2VH1 = HlV^ = HlV— = -HlV^ (23) 

til Ui ri2 

is a vector field in W^'°^{X) such that /3 ^ (on a set of positive measure), 
(iii') same as (iii) above with 

\l3\{x) > ao > 0, a.e. xeX. (24) 



10 



Beyond the regularity assumptions on (7, a, F), the domain X, and the boundary con- 
ditions /i and /2, the only real assumption we impose is (12^ . In general, there is no 
guaranty that the gradient of ^ does not vanish. In dimension d = 2, a simple condition 
guarantees that (!24|) holds. We have the following result [H HT] : 

Lemma 3.1 ( |13j ) Assume that h = ^ on dX is an almost two-to-one function in the 
sense of l4^ , i.e., a function that is a two-to-one map except possibly at its minimum 
and at its maximum. Then ([21]) is satisfied. 

In dimension > 3, the above result on the (absence of) critical points of elliptic 
solutions no longer holds. By continuity, we verify that ([21]) is satisfied for a large class 
of illuminations when 7 is close to a constant and a is sufficiently small. For arbitrary 
coefficients (7, a) in dimension c? > 3, a proof based on CGO solutions shows that (!24l) 
is satisfied for an open set of illuminations; see [18] and section 15.21 below. Note also 
that (124)) is a sufficient condition for us to solve the inverse problem of QPAT. In [1], a 
similar problem is addressed in dimension d = 2 without assuming a constraint of the 
form (121). 

We first prove a result that provides uniqueness up to a specified transformation. 

Theorem 3.2 ( |13|, 118]) Assume that hypotheses (i)-(iii) hold. Then 

(a) Hi{x) and H2{x) uniquely determine the measurement operator T-L : H^^dX) — 
H^{X), which to f defined on dX associates 'H{f) = H in X defined by 

(b) The measurement operator T-L uniquely determines the two functional: 

X(x):=^(x), g(a;):=-(^ + -)(x). (25) 
Fa ^ 7^ 

Here A is the Laplace operator. 

(c) Knowledge of the two functionals x o.iT'd q uniquely determines Hi{x) and H2{x). 
In other words, the reconstruction of (7, a, F) is unique up to transformations that 
leave {x, q) invariant. 

The proof of this theorem is given in [13] under the additional assumption (iii'). The 
following minor modification allows one to prove the theorem as stated above. The proof 
in [13] is based on the fact that JxiP ~ l)'^\(3\'^dx = implies that p = 1 a.e. Under 
assumption (iii), (3 does not vanish on a set of positive (Lebesgue) measure. However, 
P is the solution of (!26|) below with bounded from below by a positive constant. 
This implies that ^ is the solution of an elliptic equation. If /3 vanishes on a ball, and 
hence U2 = Cui for some constant C, on a set of positive measure, hence on a ball, then 
U2 — Cui is a vanishing solution of (|2T]) . and hence vanishes everywhere by the unique 
continuation principle; see e.g. [361 Chapter 3]. This would violate (iii). 

3.1.2 Reconstruction of two coefficients. 

The above result shows that the unique reconstruction of (7, a, F) is not possible even 
from knowledge of the full measurement operator "H defined in Theorem 13.21 Two well- 
chosen illuminations uniquely determine the functionals (x, q) and acquiring additional 
measurements does not provide any new information. However, we can prove that if one 
coefficient in (7, a, F) is known, then the other two coefficients are uniquely determined: 
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Corollary 3.3 ( |13] ) Under the hypotheses of the previous theorem, let {x,(l) in 
he known. Then: 

(a) IfT is known, then (7,0") are uniquely determined. 

(b) If --f is known, then (a, F) are uniquely determined, 
(a) If a is known, then (7, F) are uniquely determined. 

The above uniqueness results are constructive. In all cases, we need to solve the 
following transport equation for x'- 

- V ■ = in X, x\dx known on dX, (26) 

with /3 the vector field defined in fl23|) . This uniquely defines x > 0- Then we find that 



This provides exphcit reconstructions for (x,?) from knowledge of [Hi,H2) when (l24l) 
holds. 

In case (b), no further equation needs to be solved. In cases (a) and (c), we need to 
solve an elliptic equation for ^^7, which is the linear equation 

(A + g)V^+i- = 0, X, y7|ax = VW' -^X, (28) 
in (a) and the (uniquely solvable) nonlinear (semi-linear) equation 

y^(A + g)y^ + a = X, ^/l\ax = ^lm^ (29) 

in (c). These inversion formulas were implemented numerically in [13]. Moreover, 
reconstructions are known to be Holder or Lipschitz stable depending on the metric 
used in the stability estimate. For instance, we have: 

Theorem 3.4 ( |13] ) Assume that the hypotheses of Theorem \3.2\ and (iii') hold. Let 

H = {Hi.H-i) he the measurements corresponding to the coefficients (7, cr, F) for which 
hypothesis (iii) holds. Let H = [Hi, H2) he the measurements corresponding to the same 
illuminations (/i,/2) with another set of coefficients (7, a, F) such that (i) and (ii) still 
hold. Then we find that 

\\X - xWl^x) <C\\H- ^11^^^ for all2<p<oo. (30) 

Let us assume, moreover, that 7(x) is of class C^{X). Then we have that 

\\X - xlUoc(x) < C\\H - H\\f^ for all2<p< 00. (31) 

We may for instance choose p = 4 above to measure the noise level in the measurement 
H in the square integrable norm when noise is described by its power spectrum in the 
Fourier domain. This shows that reconstructions in QPAT are Holder stable, unlike the 
corresponding reconstructions in Optical Tomography fLOl [60] . 
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3.1.3 An application to Transient Elastography 



We can apply the above results to the time-harmonic reconstruction in a simplified model 
of transient elastography. Let us assume that 7 and p are unknown functions of x G X 
and a; G M. Recall that the displacement solves ( TTH]) . Assuming that u{x]uj) is known 
after step 1 of the reconstruction using the ultrasound measurements, then we are in 
the setting of Theorem 13 . 2 1 wit h Ta = 1. Let us then assume that the two illuminations 
fi{x;u) and f2{x]u) are chosen such that for ui and U2 the corresponding solutions of 
( 1T8|) . we have that holds. We have seen a sufficient condition for this to hold in 
dimension n = 2 in Lemma 13.11 and will present other sufficient conditions in section 
15.21 below devoted to CGO solutions in the setting n > 3. Then, fl2^ shows that the 
reconstructed function x uniquely determines the Lame parameter 'j{x; u) and that the 
reconstructed function q then uniquely determines u'^p and hence the density parameter 
p{x;uj). The reconstructions are performed for each frequency uj independently. We 
may summarize this as follows: 

Corollary 3.5 Under the hypotheses Theore'm \3.2\ and the hypotheses described above, 
let (XtQ) in (!25|) be known. Then {'~f{x;uj), p{x;uj)) are uniquely determined by two 



well-chosen measurements. Moreover, the stability results in Theorem\3.4\ hold. 



Alternatively, we may assume that in a given range of frequencies, 7(x) and p{x) 
are independent of u. In such a setting, we expect that one measurement u{x; u) for 
two different frequencies will provide sufficient information to reconstruct {'y{x), p{x)). 
Assume that u{x; u) is known for u = Uj, j = 1,2 and define < a = oj\uj^'^ 7^ 1. Then 
straightforward calculations show that 

V-7/3„ = 0, /3« = (uiVu2-au2Vui). (32) 

This provides a transport equation for 7 that can be solved stably provided that > 
Co > 0, i.e., /3o does not vanish on X. Then, Theorem 13.21 and Theorem 13.41 apply in 
this setting. Since cannot be written as the ratio of two solutions as in fl23|) when 
Q! = 1, the results obtained in Lemma [3. II do not apply when a ^ 1. However, we prove 
in section that ^ Cq > is satisfied for an open set of illuminations constructed 
by means of CGO solutions for all a > 0; see (!96|) below. 



3.1.4 Reconstruction of one coefficient 

Let us conclude this section by some comments on the reconstruction of a single co- 
efficient from a measurement linear in u. From an algorithmic point of view, such 
reconstructions are siginficantly simpler. Let us consider the framework of Corollary 
13.31 When F is the only unknown coefficient, then we solve for u in fl2ip and reconstruct 
F from knowledge of H. 

When only cr is unknown, then we solve the elliptic equation for u 

—V ■ 7VM + — = in A, u = g on dX, 

and then evaluate a = 

i u 
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When only 7 is unknown with either H = au in QPAT or with H = u'm elastography 
or in apphcations to ground water flows [HET], then u is known and 7 solves the following 
transport equation 

—V ■ 7Vn = 5* in X, 7 = 7|gx on 9X, 

with 5* known. Provided that the vector field Vm does not vanish, the above equation 
admits a unique solution as in fl26|) . The stability results of Theorem 13.41 then apply. 
Other stability results based on solving the transport equation by the method of char- 
acteristics are presented in [51]. In two dimensions of space, the constraint that the 
vector field Vm does not vanish can be partially removed. Under appropriate conditions 
on the oscillations of the illumination g on dX, stability results are obtained in [1] in 
cases where Vm is allowed to vanish. 

3.2 Reconstructions from quadratic functionals in u 
3.2.1 Reconstructions under smallness conditions 

The TAT and (simplified) AOT problems are examples of a more general class we define 
as follows. Let D) be an operator acting on functions defined in for m e N* 
an integer and with values in the same space. Consider the equation 

P{x,D)u = a{x)u, xeX 
u = f, X e dX 

We assume that the above equation admits a unique weak solution in some Hilbert space 
T-Li for sufficiently smooth illuminations f{x) on dX. 

For instance, P could be the Helmholtz operator ik~^{A + k"^) seen in the preceding 
section with u G Tii := H^{X;C) and / G H2{dX;C). Time-harmonic Maxwell's 
equations can be put in that framework with m = n and 

P(x, D) = i-(V X V X -e). (34) 
ik 

We impose an additional constraint on P{x, D) that the equation P(x, D)u = / on 
X with u = on dX admits a unique solution m % = L^(X;C'"). For instance, 
l-L = L'^{X] C) in the example seen in the preceding section in the scalar approximation 
provided that k"^ is not an eigenvalue of —A on X. For Maxwell's equations, the above 
constraint is satisfied so long as k'^ is not an internal eigenvalue of the Maxwell operator 
[21]. This is expressed by the existence of a constant a > such that: 

(P(x, D)u, u)h > «(m, u)h- (35) 

We assume that the conductivity a is bounded from above by a positive constant: 

< a{x) < aM a-e. x G X. (36) 

We denote by Em the space of functions a{x) such that fl36|) holds. Measurements are 
of the form H{x) = a{x)\u\'^, where | ■ | is the Euclidean norm on C™. Then we have 
the following result. 
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(33) 



Theorem 3.6 ( |16] ) Let aj G forj = 1,2. Let Uj be the solution to P{x,D)uj = 
(JjUj in X with uj = f on dX for j = 1,2. Define the internal Junctionals Hj{x) = 
aj{x)\uj{x)\'^ on X . 

Then for ctm sufficiently small so that < a, we find that: 
(i) [Uniqueness] If Hi = H2 a.e. in X , then ai{x) = (x) a.e. in X where Hi = H2> 0. 
(a) [Stability] Moreover, we have the following stability estimate 

\\{V^l-^/^2)Wl\\n < C\\{y^i-y^2)w2\\H, (37) 
for some universal constant C and for positive weights given by 

Wt(x) = TT -^^^(x), W2(x) = ^ — , max ^^/—^ (x) + max —=r- (x) . (38) 

•ii2^i «-suPv/aIa^i=i:2 ' j=i,2 ^ 

Here j' = is defined as = 2 and j'(2) = 1. 

The theorem uses the spectral gap in ( !35l) . Some straightforward algebra shows that 



P{x,D){ui - U2) = ^/o^{\u2\ul - \u1\u2) + (v^- V^)(-r^ 



P2 



Here we have defined it = Although this does not constitute an equation for U1 — U2, 
it turns out that 

This combined with ( l35l) yields the theorem after some elementary manipulations [T6] . 



3.2.2 Reconstructions for the Helmholtz equation 

Let us consider the scalar model of TAT. We assume that a G H^{X) for p > ^ and 
construct 

q(x) = e + ika{x) G H^iX), p > -. (39) 

We assume that q{x) is the restriction to X of the compactly supported function (still 
called q) q E H'^iW^). The extension is chosen so that |T8] ||g|x||//p(x) < C'||g||i^p(iRn) for 
some constant C independent of q. Then (|6]) may be recast as 

Am + q{x)u = in X, u = f on dX. (40) 

The measurements are of the form H[x) = a{x)\u\'^{x). 

The inverse problem consists of reconstructing q{x) from knowledge of H{x). Note 
that q{x) need not be of the form f l59]) . It could be a real- valued potential in a Helmholtz 
equation as considered in [59] with applications in the so-called inverse medium problem. 
The reconstruction of q{x) in ( 140|) from knowledge of H{x) = (j{x)\u\^{x) has been 
analyzed in [161 EH] • 

In the first reference, we have the following global stability reconstruction result. We 
define Y = H'p{X) for p > f and M. as the space of functions in Y with norm bounded 
by a fixed (arbitrary) M > 0. Let us define Z = H^'i (dX). 
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Theorem 3.7 ( [16] ) Let a and a he functions in M.. Let f E Z be a given (complex- 
valued) illumination and H{x) be the measurement given in ([7j) for u solution of (|6]). 
Let H{x) be the measurement constructed by replacing a by a in and IQ. 

Then there is an open set of illuminations f in Z such that H{x) = H{x) in Y 
implies that a{x) = (t{x) in Y. Moreover, there exists a constant C independent of a 
and a in M. such that 

\\a -a\\Y <C\\H - H\\y. (41) 

The theorem is written in terms of a, which is the parameter of interest in TAT. The 
same result holds if a is replaced by q{x) in (141 p . The reconstruction of a is also 
constructive as the application of a Banach fixed point theorem. The proof is based on 
the construction of complex geometric optics solutions that will be presented in section 

m 

In reference [SH], the following local stability result is obtained. 

Theorem 3.8 ( |59] ) Let q{x) > cq > be real-valued, positive, bounded on X and 
such that is not an eigenvalue o/ A + g with domain Hq{X) fl H2{X). Let q satisfy 
the same hypotheses and let H and H be the corresponding measurements. 

Then there is a constant e > such that if q and q are e— close in L°°(X) and if f is 
in an e— dependent open set of (complex-valued) illuminations, then there is a constant 
C such that 

\\q-q\\mx)<C\\H-H\\L2^x)- (42) 

Both theorems 13.71 and 13.81 show that the TAT and the inverse medium problem are 
stable inverse problems. This is confirmed by the numerical reconstructions in [16]. The 
first result is more global but requires more regularity of the coefficients. It is based 
on the use of complex geometric optics CGO solutions to show that an appropriate 
functional is contracting in the space of continuous functions. The second result is 
more local in nature (a global uniqueness result is also proved in [59]) but requires less 
smoothness on the coefficient q{x) and provides a stability estimate in the larger space 
L'^{X). It also uses CGO solutions to show that the norm of a complex- valued solution 
to an elliptic equation is bounded from below by a positive constant. In both cases, the 
CGO solutions have traces at the boundary dX and the chosen illumination / needs to 
be chosen in the vicinity of such traces. 

The results obtained in Theorem 13.61 under smallness constraints on a apply for 
very general illuminations /. The above two results apply for more general (essentially 
arbitrary) coefficients but require more severe constraints on the illuminations /. 

3.2.3 Non-unique reconstruction in the AOT setting 

The above results concern the uniqueness of the reconstruction of the potential in a 
Helmholtz equation when well-chosen complex-valued boundary conditions are imposed. 
They also show that the reconstruction of < cq < q{x) in Au -\- qu = with real- 
valued u = f from knowledge of qu"^ is unique. This corresponds to P{x, D) = —A. 
In a simplified version of the acousto-optics problem considered in [T7], it is interesting 
to look at the problem where P{x, D) = A and where the measurements are given by 
H{x) = (j(x)m^(x). Here, u is the solution of the elliptic equation (—A + a)u = on X 
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with u = / on dX. Assuming that / is non-negative, which is the physically interesting 
case, we obtain that \u\ = u and hence 



A(ni - U2) = y/o\o^{u2 - ui) + - V^) (- 



CTi (72 



'Hi VH2 

Therefore, as soon as is not an eigenvalue of A + s/aid^, we obtain that Ui = U2 and 
hence that ai = 02- For ctq such that is not an eigenvalue of A + (Tq, we find that for 
a\ and 02 sufficiently close to ctq, then = -f^2 implies that (j\ = (J2 on the support of 
Hi = H2. 

However, it is shown in [15] that two different, positive, absorptions ctj for j = 1, 2, 
may in some cases provide the same measurement H = ajuj with Auj = ajUj on X 
with Uj = / on dX and in fact cxi = a2 on dX so that these absorptions cannot be 
distinguished by their traces on dX. This counter-example shows that conditions such 
as the smallness condition in Theorem 13.61 are necessary in general. 

More generally, and following [T3], consider an elliptic problem of the form 

Pu = au in X, u = f on dX, (43) 

and assume that measurements of the form H{x) = a{x)u'^{x) are available. Here, 
P is a self-adjoint, non-positive, elliptic operator, which for concreteness we will take 
of the form Pu = V ■ 7(x)Vm with 7(x) known, sufficiently smooth, and bounded 
above and below by positive constants. We assume / > and 7 > so that by the 
maximum principle, m > on X. We also assume enough regularity on dX and / so 
that u e C^'f^iX) for some /3 > [SU]. 
We observe that 

uPu = H in X, u = f on dX (44) 

so that the inverse problem may be recast as a semilinear problem. The non-uniqueness 
result is an example of an Ambrosetti-Prodi result [3] and in some sense generalizes the 
observation that x — admits 0, 1, or 2 (real- valued) solution(s) depending on the 
value of x^. Let us define 

: C^'^iX) C°'^{X), u H> 0(m) = uPu. (45) 

The singular points of are calculated from its first-order Frechet derivative: 

(j)'{u)v = vPu + uPv. (46) 

The operator 0'(m) is not invertible when cr := ^ is such that P + Xa admits A = 1 
as an eigenvalue. Let ctq be such that P + ao is not invertible. We assume that the 
corresponding eigen-space is one dimensional and spanned by the eigenvector > on 
dX such that (P + ctq)^ = and = on dX. Let us define uq as 

Puq = (JqUq in X, Uq = f on dX, ctq > 0. (47) 

Moreover, uq is a singular point of (f){u) with (j)'{uQ)ip = . Then define 

us := Uq + 6ij, X, 6 e (-5o, So) (48) 

as := = (To Hs := asug = aousu^s = cro{uQ- d f ). (49) 

Us uo + dip 

We choose 5o such that cr^ > a.e. on X for all 5 G (— 5o, Then, we have: 
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Proposition 3.9 ( |15] ) Let Uq be a singular point and Hq = a critical value of 

(p as above and let ip be the normalized solution of (j)' {uo)'iIj = 0. Let us, as, and Hs be 
defined as in (I48|) - (l49|) for 7^ 5 G (— ^o^^o) for 5q sufficiently small. Then we verify 
that: 

as 7^ as > 0, Hs = H^s, Pus = asus in X, us = f on dX. 

This shows the non-uniqueness of the reconstruction of a from knowledge oi H = au^. 
Moreover we verify that a±s agree on dX so that this boundary information cannot 
be used to distinguish between as and a^s- The non-uniqueness resuh is not very 
restrictive since we have seen that two coefficients, hence one coefficient, may be uniquely 
reconstructed from two well-chosen illuminations in the PAT results. Nonetheless, the 
above result shows once more that identifiability of the unknown coefficients is not 
always guaranteed by the availability of internal measurements. 



4 Reconstructions from functionals of Vu 

We have seen two models of hybrid inverse problems with measurements involving Vu. 
In UMEIT, the measurements are of the form H{x) = 7(x)|Vup(x) whereas in CDII, 
they are of the form H{x) = '-f{x)\'Vu\{x). 

Let us consider more generally measurements of the form H{x) = 7(a;)|VMp~^ for u 
the solution to the elliptic equation 

- V ■ 7(x)Vm = in X, u = f on dX. (50) 

Since H{x) is linear in 7(x), we have formally what appears to be an extension to p > 
of the p— Laplacian elliptic equations 

-V-^^Vn = 0, (51) 

posed on a bounded, smooth, open domain X C M", n > 2, with prescribed Dirichlet 
conditions, say. When 1 < p < 00, the above problem is known to admit a variational 
formulation with convex functional J[Vm] = H{x)\Vu\Pdx, which admits a unique 
minimizer in an appropriate functional setting solution of the above associated Euler- 
Lagrange equation [26] . 

When p = 1, the equation becomes degenerate while for p < 1, the equation is 
in fact hyperbolic. We consider the problem of measurements that are quadratic (or 
bilinear) in Vn (with applications to UMEIT and UMOT) in the next two sections. In 
the following section, we consider the case p = 1. 



4.1 Reconstruction from a single power density measurement 

The presentation follows [1]. When p = so that measurements are of the form H{x) = 
7(x)|VMp, the above 0— Laplacian turns out to be a hyperbolic equation. Anticipating 
this behavior, we assume the availability of Cauchy data (i.e., u and 7Z/ ■ Vu with u the 
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unit outward normal to X) on dX rather than simply Dirichlet data. Then f lSip with 
p = becomes after some algebra 

(/ - 2Vn ® Vu) : V\ + V In /J ■ Vn = in X, u = f and -7^ = J on dX. (52) 



du 



Here Vm = With 



^^^^ = ^^^(Vu) = -S'^ + 2(Vu)i(Vn)j and F = -(VlniJ)^, (53) 
the above equation is in coordinates 

g'\Vu)d^jU + k'diU = inX, u = f and ^=J on ^X. (54) 

Since g^^ is a definite matrix of signature {l,n — 1), then fl5^ is a quasilinear strictly 
hyperbolic equation. The Cauchy data / and j then need to be provided on a space-like 
hyper-surface in order for the hyperbolic problem to be well-posed [33] . This is the main 
difficulty with solving (l54l) with redundant Cauchy boundary conditions. 

In general, we cannot hope to reconstruct u{x), and hence 7(x) on the whole domain 
X. The reason is that the direction of "time" in the second-order hyperbolic equation 
is Vm(x). The normal z/(x) at the boundary dX will distinguish between the (good) 
part of dX that is "space-like" and the (bad) part of dX that is "time-like" . Space-like 
surfaces such as t = provide stable information to solve the standard wave equation 
whereas in general it is known that arbitrary singularities can form in a wave equation 
from information on "time-like" surfaces such as a; = or y = in a three dimensional 
setting (where {t,x,y) are local coordinates of X) [33] . 

4.1.1 Uniqueness and stability 

Let (m, 7) and {u, 7) be two solutions of the Cauchy problem (IMI) with measurements 
[H,f,j) and {H,f,j), where we define the reconstructed conductivities 

H H 

Let V = u — u. We find that 

V-( |^,iff^ i„ |(VM+V72)0(V^+V72)-(|V^|'+|V72nj|Vt;+^i/(-^+-^)) = 0. 
V|Vm| |Vm| I- J V|Vm| |Vm| // 

This equation is recast as 



Q'^{x)dlv + rd,v + d,{r5H) = Q inX, v = f-f, ^=3-3 on dX, (56) 



for appropriate coefficients and where 

^ (Vm + Vm) ® (Vm + Vm) - (|Vm|^ + \Vu\^)I } 

(57) 



g(x) = rr=^|(VM + Vm) ® (Vm + Vm) - (|VMp + |VMn/| 
= a{x) ( e(x) ® e(x) — /3^(x) (/ — e(x) ® e(x)) 
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with 

, , Vn + Vu 2 |Vn + Vnp 

|Vn + Vu| I Vu + Vn|^ — (I Vnp + I Vupj 

and a{x) is the appropriate (scalar) normahzation coefficient. For Vu and V£t suf- 
ficiently close so that Vn ■ Vn > 0, then the above linear equation for v is strictly 
hyperbolic. We define the Lorentzian metric {) = so that are the coordinates 
of the inverse of the matrix q^K We denote by (■, ■) the bilinear product associated to 
f) so that (u, v) = \)iju'^v^ where the two vectors u and v have coordinates and f \ 
respectively. We verify that 

The space-like part Eg of dX is given by {)(z/, v) > 0, i.e., i/ is a time-like vector, or 
equivalently 

IK^)-e(x)|2>^^^ xedX. (60) 

Above, the "dot" product is with respect to the standard Euclidean metric and z/ is a 
unit vector for the Euclidean metric, not for the metric f). Let Si be an open connected 
component of and let O = Uo<r<sS2(r) be a domain of infiuence of Ei swept out 
by the space- like surfaces S2(r); see [9l[58]. Then we have the following local stability 
result: 

Theorem 4.1 (Local Uniqueness and Stability.) Let u and u he two solutions of 
( 15^ . We assume that q constructed in ( 1571) is strictly hyperbolic. Let Ei he an open 
connected component of E^ the space-like component of dX and let O he a domain of 
infiuence of Ei constructed as ahove. Let us define the energy 

E{dv) = {dv,U2)^ - ^{dv,dv){l^2,l^2). (61) 

Here, dv is the gradient of v in the metric \) given in coordinates hy g^^djv. Then: 

[ E{dv)dx<c( I + f \\/5H\'^dx), (62) 

Jo ^ Jsi Jo ' 

where dx and da are the standard measures on O and Ei, respectively. 

In the Euclidean metric, let V2{x) he the unit vector to x E E2(t), define c(x) : = 
1^2 (x) ■ e(x) and 

(63) 



6 := min 



cHx) - ^ 



l+/32(x)J 



Then we have that 



I \v'\ + \Vv\' + {j-^fdx<^n \f-fy + \j-j\^dj+ f \V6H\'dx), (64) 

where 7 and 7 are the conductivities in fl55|) . Provided that f = f, j = j, and H = H , 
we ohtain that u = and the uniqueness result u = u and 7 = 7. 
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The proof is based on adapting energy methods for hyperbohc equations as they are 
summarized in [58]. The energy E{dv) fails to control dv for null-like or space-like 
vectors, i.e., i){dv, dv) < 0. The parameter 6 measures how time-like the vector dv is on 
the domain of influence O. As O approaches the boundary of the domain of influence 
of T,g and 9 tends to 0, the energy estimates deteriorate as indicated in 

Assuming that the errors on the Cauchy data / and j are negligible, we obtain the 
following stability estimate for the conductivity 

C 

h - ih^o) < "^11-^ ~ H\\hi(x)- (65) 

Under additional regularity assumptions on 7, for instance assuming that H G H^[X) 
for s > 2, we find by standard interpolation that 

||7 - iWlho) <^\\H- H\\%)\\H + ^lli(x), (66) 

We thus obtain Holder-stable reconstructions in the practical setting of square integrable 
measurement errors. However, stability is local. Only on the domain of influence of the 
space-like part of the boundary can we obtain a stable reconstruction. This can be 
done by solving a nonlinear strictly hyperbolic equation analyzed in [H] using techniques 
summarized in [M] . 

4.1.2 Global reconstructions 

In the preceding result, the main roadblock to global reconstructions was that the do- 
main of influence of the space-like part of the boundary was a strict subset of X . There 
is a simple solution to this problem: simply make sure that the whole boundary is a 
level set of u and that no critical points of u (where Vu = 0) exist. Then all of X is in 
the domain of influence of the space-like part of 9X, which is the whole of dX. This 
setting can be made possible independent of the conductivity 7 in two dimensions of 
space but not always in higher dimensions. 

Let n = 2. We assume that X is an open smooth domain diffeomorphic to an annulus 
with boundary dX = OXq U dXi . We assume that / = on the external boundary OXq 
and that / = 1 on the internal boundary dXi. The boundary of X is thus composed of 
two smooth connected components that are different level sets of the solution u. The 
solution u to fl50|l is uniquely defined on X. Then we can show: 

Proposition 4.2 ([9]) We assume that both the geometry of X and'j{x) are sufficiently 
smooth. Then |Vu| is bounded from above and below by positive constants. The level 
sets Sc = {x G X, u{x) = c} for < c < 1 are smooth curves that separate X into two 
disjoint subdomains. 

The proof is based on the fact that critical points of solutions to elliptic equations in two 
dimensions are isolated [1]. The result extends to higher dimensions provided that |Vn| 
does not vanish with exactly the same proof. In the absence of critical points, we thus 
obtain that e(x) = Vn = z/(x) so that z/(x) is clearly a time-like vector. Then the local 
results of Theorem 14.11 become global results, which yields the following proposition: 
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Proposition 4.3 Let X be the geometry described above in dimension n>2 and u{x) 
the solution to (!50|l . We assume here that both the geometry and ■j{x) are sufficiently 
smooth. We also assume that \Wu\ is bounded from above and below by positive constants. 
Then the nonlinear equation (154|) admits a unique solution and the reconstruction of u 



and of 'J is stable in X in the sense described in Theorem\4.1 



In dimensions n > 3, we cannot guaranty tliat u does not liave any critical point 
independent of the conductivity. If the conductivity is close to a constant, then by 
continuity, u does not have any critical point and the above result applies. This proves 
the result for sufficiently small perturbations of the case 7(x) = 7o- In the general 
case, however, we cannot guaranty that Vn does not vanish and in fact can produce 
counter-examples (see [9]): 

Proposition 4.4 ([9|, I20|, 145] ) There is an example of a smooth conductivity such that 
u admits critical points. 

So in dimensions n > 3, we are not guaranteed that the nonlinear equation will remain 
strictly hyperbolic. What we can do, however, is again to use the notion of complex 
geometric optics solutions. We have the result: 

Theorem 4.5 ([9]) Let 7 be extended by-fo = 1 onM"\X, where X is the domain where 
7 is not known. We assume that 7 is smooth on M". Let 7(x) — 1 be supported without 
loss of generality on the cube (0, l)x(— i, ^)"~^. Define the domainX = (0, l)xBn-i{a), 
where 5„_i(a) is the n — 1-dimensional ball of radius a centered at and where a is 
sufficiently large that the light cone for the Euclidean metric emerging from Bn-i{a) 
strictly includes X. Then there is an open set of illuminations (/i,/2) such that if ui 
and U2 are the corresponding solutions of (l50l) . then the following measurements 

=7(x)|VmiP(x), H2{x)=^{x)\Vu2\\x), H:,{x) = ^{x)\V{u^ + U2)\\ (67) 

with the corresponding Cauchy data (/i,ji), (/2,i2) (ind (/i + /2,ji + ^2) at Xi = 
uniquely determine 7(x). Moreover, let Hi be measurements corresponding to 7 and 
(/15 Ji) anc? (/2, J2) the corresponding Cauchy data at xi = 0. We assume that 7(x) — 1 
and 7(0;) — 1 are smooth and such that their norm in 2+^+^(M'") for some e > are 
bounded by M. Then for a constant C that depends on M , we have the global stability 
result 

3 

||7 - 7lL2(x) < C(\\dc - rfc||(L2(B„.i(a)))4 + 11^^* - Vi/.||L2(X)). (68) 

1=1 

Here, we have defined dc = (/i, Jij /27 J2) with dc being defined similarly. 

The "three" measurements Hi in (l67j) actually correspond to "two" physical measure- 
ments since H^ may be determined from the experiments yielding Hi and H2 as demon- 
strated in P E]. The three measurements are constructed so that two independent 
strictly hyperbolic Lorentzian metrics can be construct everywhere inside the domain. 
These metrics are constructed by means of CGO solutions. The boundary conditions fj 
have to be close to the traces of the CGO solutions. We thus obtain a global Lipschitz 
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stability result. The price to pay is that the open set of illuminations is not very explicit 
and may depend on the conductivities one seeks to reconstruct. 

For conductivities that are close to a constant, several reconstructions are therefore 
available. We have seen that geometries of the form of an annulus (with a hole that 
can be arbitrarily small and arbitrarily close to the boundary where / = 0) allowed us 
to obtain globally stable reconstructions since in such situations, it is relatively easy 
to avoid the presence of critical points. The method of CGO solutions can be shown 
to apply for a well-defined set of illuminations since the (harmonic) CGO solutions are 
explicitly known for the Euclidean metric and of the form e'''^ for p a complex valued 
vector such that p ■ p = 0. After linearization in the vicinity of the Euclidean metric, 
another explicit reconstruction procedure was introduced in [H]. 

4.2 Reconstructions from multiple power density measurements 

Rather than reconstructing 7 from one given measurement of the form 7(x)|VMp, we 
can instead acquire several measurements of the form 

Hij{x) = 7(x)VMi(x) ■ Vuj{x) in X, 1 < i, j < M, (69) 

where Uj solves the elliptic problem f lSU]) with / given by fj for 1 < j < M. The result 
presented in Theorem 14 . 5 1 above provides a positive answer for M = 2 when the available 
internal functionals are augmented by Cauchy data at the boundary of the domain of 
interest. 

Results obtained in [TTl [211 HSj and based on an entirely different procedure and not 
requiring knowledge of boundary data show that M = 2[^^J measurements allow for a 
global reconstruction of 7, i.e., M = n for n even and M = n + 1 for n odd. Such results 
were first obtained in [21] in the case n = 2 and have been extended with a slightly 
different presentation to the cases n = 2 and = 3 in [11] while the general case n > 2 
is treated in [16]. Let us assume that n = 3 for concreteness. Then Hij = Si ■ Sj, where 

we have defined 

Sj{x) = ^/-f{x)Wuj{x), 1<J<M. 

Let S = {Si, S2, S3) be a matrix of n = 3 column vectors Sj. Then S'^S = H, where 
S'^ is the transpose matrix made of the rows given by the Sj. We do not know S or the 
Sj, but we know its normal matrix S^S = H. Let T be a matrix such that R = ST'^ 
is a rotation- valued field on X. Two examples are T = or the lower-triangular T 
obtained by the Gram-Schmidt procedure. We thus have information on S. We need 
additional equations to solve for S, or equivalently R, uniquely. The elliptic equation 
may be written as V ■ y/ySj = 0, or equivalently 

V-S, + F-S, = 0, F = V(log v/^) = ^Vlog7. (70) 

Now, since 7~^>S'j is a gradient, its curl vanishes and we find that 

V X Sj - F X Sj = 0. (71) 

Here, F is unknown. We first eliminate it from the equations and then find a closed 
form equation for S or equivalently for i? as a field in SO{n; M). 
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Let T be the aforementioned matrix T, say T = H 2 with entries tij for 1 < i, j < n. 
Let be the entries of and define the vector fields 

Vij := V{Uk)t''', i.e., l^A = 1 < ^, J, / < n. (72) 

We then define R{x) = S{x)T'^{x) G SO{n; M) the matrix whose columns are composed 
of the column vectors Rj = SjT^. Then in all dimension n > 2, we find 

Lemma 4.6 146] ) In n > 2, we have the following expression: 

F = -(-Vlogdeti/ + J2 ((^^^- + ^^•*) ■ (73) 

The proof in dimension n = 2,3 can be found in [11] and in arbitrary dimension in [l6]. 

Note that the determinant of H needs to be positive on the domain X in order for 
the above expression for F to make sense. It is, however, difficult to ensure that the 
determinant of several gradients remains positive and there are in fact counter-examples 
as shown in [20]. Here again, complex geometric optics solutions are useful to control the 
determinant of gradients of elliptic solutions locally and globally using several solutions. 
We state a global result in the practical setting n = 3. 

Let be m > 3 solutions of the elliptic equation and assume that there exists an open 
covering O = {fik}i<k<N C U^^fi.;), a constant Cq > and a function r : [1, A^] 3 
i H- r(i) = (r(i)i, r(i)2, r(i)3) G [l,m]^, such that 

inf det{Sr{i),{x), 5^(^)2(0;), 5^(^)3(0;)) > cq, l<i< N. (74) 
Then we have the following result: 

Theorem 4.7 (3D global uniqueness and stability) Let X G M.^ be an open con- 
vex bounded set, and let two sets of m > 3 solutions of (!50|) generate measurements 
{H,H) whose components belong to W^'^{X). Assume that one can define a couple 
{O, t) such that ( 171|) is satisfied for both sets of solutions S and S . Let also Xq G VLi^ C X 
and'j{xQ), 7(xo), {ST{io)i{xo), Sr{io)i{^o)}i<i<3 be given. Let'y and'j be the conductivities 
corresponding to the measurements H and H , respectively. Then we have the following 
stability estimate: 

II log7 - log7||H/i.-(x) < C{eQ + lli^" - ^||h/i,-(x)), (75) 

where eo is the error at the initial point xq 

3 

eo = I log 70 - log7o| + ^ \\ST{io)Axo) - Sr{io)Axo)\\- 

1=1 

This shows that the reconstruction of 7 is stable from such redundant measurements. 
Moreover, the reconstruction is constructive. Indeed, after eliminating F from the 
equations for R, we find an equation of the form Vi? = G{x,R), where G{x,R) is 
polynomial of degree three in the entries of R. This is a redundant equation whose 
solution, when it exists, is unique and stable with respect to perturbations in G and the 
conditions at a given point xq. 

That fTMl) is satisfied can again be proved by means of complex geometric optics 
solutions as is briefiy mentioned in section [5l2] below: see [TT] . 
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4.3 Reconstruction from a single current density measurement 

Let us now come back to the 1— Laplacian, which is a degenerate elhptic problem. In 
many cases, this problem admits multiple admissible solutions [38] . The inverse problem 
then cannot be solved uniquely. In some settings, however, uniqueness can be restored 

[Mmzmsiiis]. 

Recall that the measurements are of the form H[x) = 7|Vm| so that u solves the 
following degenerate quasi-linear equation 

V • ^^Vn = in X. (76) 
I Vm| 

Different boundary conditions may then be considered. It is shown in that the 
above equation augmented with Neumann boundary conditions of the form 

H Ou 

h on dX, I uda = 0, 



\Vu\diy Jqx 

admits an infinite number of solutions once it admits a solution, and may also admit 
no solution at all. One possible strategy is to acquire two measurements of the form 
H{x) = 7|Vu| corresponding to two prescribed currents. In this setting, it is shown 
in that (appropriately defined) singularities of 7 are uniquely determined by the 
measurements. We refer the reader to the latter reference for the details. 

Alternatively, we may augment the above equation (1761) with Dirichlet data. Then 
the reconstruction of 7 was shown to be uniquely determined in [ITJ HHl HQ]. Why 
Dirichlet conditions help to stabilize the equation may be explained as follows. The 
1— Laplace equation (|75]) may be recast as 

(/ - Vu O Vu) : V^n + V In if ■ Vn = 0, 



following similar calculations to those leading to fl52|) . The only difference is the "2" in 
front of Vu ® Vu replaced by "1" , or more generally 2 — p for a p— Laplacian. When 
p > 1, the problem remains strictly elliptic. When p < 1, the problem is hyperbolic, 
and when p = 1, it is degenerate in the direction Vu and elliptic in the transverse 
directions. We can therefore modify u so that its level sets remain unchanged and still 
satisfy the above partial differential equation. This modification can also be performed 
so that Neumann boundary conditions are not changed. This is the procedure used in 
|38] to show the non-uniqueness of the reconstruction for the 1-Laplacian with Neumann 
boundary conditions. 

Dirichlet conditions, however, are modified by changes in the level sets of u. It turns 
out that even with Dirichlet conditions, several (viscosity) solutions to fl76|) may be 
constructed when H = 1; see, e.g., [ITlllS]. However, such solutions involve vanishing 
gradients on sets of positive measure. 

The right formulation for the CDII inverse problem that allows one to avoid vanishing 
gradients is to recast (1761) as the minimization of the functional 

F[\/v] = I H{x)\Vv\dx, (77) 
Jx 
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over V G H^{X) with v = f on dX. Note that F[Vf] is convex although it is not 
strictly convex. Moreover, let 7 be the conductivity and H = 7|Vm| the corresponding 
measurement. Let then v G H^{X) with v = f on dX. Then, 

F[Vv]= I -i\Vu\\Vv\dx> I -iVu-Vvdx= I a^fds = F[Vu], 
Jx Jx JdX 

by standard integrations by parts. This shows that u minimizes F. We have the 
following result: 

Theorem 4.8 ([49]) Let {f,H) e C^^'^idX) x C"(X) with H = 7|Vm| for some 7 G 
C°(X). Assume that H{x) > a.e. in X . Then the minimization of 

argmm {F[Vvl v G W^^\X) n C{X), , v\9x = /}, (78) 

has a unique solution uq. Moreover ao = H\Vuq\~^ is the unique conductivity associated 
to the measurement H{x). 

It is known in two dimensions of space that H{x) > is satisfied for a large class of 
boundary conditions f{x); see Lemma [5.21 in the next section. In three dimensions of 
space, however, critical points of u may arise as observed earlier in this paper; see e.g., 
[9]. The CGO solutions that are analyzed in the following section allow us to show that 
H{x) > holds for an open set of illuminations / at the boundary of the domain dX; 
see fl95l) below. 

Several reconstruction algorithms have been devised in [38| W7\ HH HQ] , to which we 
refer for additional details. The numerical simulations presented in these papers show 
that when uniqueness is guaranteed, then the reconstructions are very high resolution 
and quite robust with respect to noise in the data, as is expected for general hybrid 
inverse problems. 



5 Qualitative properties of forward solutions 
5.1 The case of two spatial dimensions 

Several explicit reconstructions obtained in hybrid inverse problems require that the 
solutions to the considered elliptic equations satisfy specific qualitative properties such 
as the absence of any critical point or the positivity of the determinant of gradients of 
solutions. Such results can be proved in great generality in dimension n = 2 but do not 
always hold in dimension n > 3. 

In dimension n = 2, the critical points of u (points x where Vm(x) = 0) are nec- 
essarily isolated as is shown in, e.g., [1]. From this and techniques of quasiconformal 
mappings that are also restricted to two dimensions of space, we can show the following 
results. 

Lemma 5.1 ([2]) Let Ui and U2 he the solutions of fl50p on X simply connected with 
boundary conditions fi = xi and f2 = X2 on dX , respectively, where x = {xi,X2) are 
Cartesian coordinates on X. Assume that 7 is sufficiently smooth. Then (xi,X2) ^ 
{ui,U2) from X to its image is a diffeomorphism. In other words, det(Vui, Vn2) > 
uniformly on X. 
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This result is useful in the analysis of UMEIT and UMOT in the case of redundant 
measurements. It is shown in [20] that the appropriate extension of this result is false 
in dimension n > 3. 

We recall that a function continuous on a simple closed contour is almost two-to-one 
if it is two-to-one except possibly at its maximum and minimum |17]. Then we have, 
quite similarly to the result in Lemma [XT] and Proposition 14. 2[ which also use the results 
in [1], the following: 

Lemma 5.2 (pTT]) Let X be a simply connected planar domain and let u he solution 
of flSOp with f almost two-to-one and a sufficiently smooth. Then \Vu\ is bounded from 
below by a positive constant on X. Moreover, the level sets of u are open curves inside 
X with their two end points on dX . 

This shows that for a large class of boundary conditions with one maximum and one 
minimum, the solution u cannot have any critical point in X. On an annulus with 
boundaries equal to level sets of u, we saw in Proposition 14.21 that u had no critical 
points on X in dimension n = 2. This was used to show that the normal vector to the 
level sets of u always forms a time- like vector for the Lorentzian metric defined in flM|) . 

All these results no longer hold in dimension n > 3. See, e.g., [HI [20] for counter- 
examples. In dimension n > 3, the required qualitative properties cannot be obtained 
for a given set of illuminations (boundary conditions) independent of the conductivity. 
However, for conductivities that are bounded (with an arbitrary bound) in an appro- 
priate norm, there are open sets of illuminations that allow us to obtain the required 
quahtative properties. One way to construct such solutions is by means of the complex 
geometric optics solutions that are analyzed in the next section. 

5.2 Complex Geometric Optics solutions 
5.2.1 CGO solutions and Helmholtz equations. 

Complex geometrical optics (CGO) solutions allow us to treat the potential q in the 
equation 

(A + g)n = in X, u = f on dX, (79) 

as a perturbation of the leading operator A. When q = 0, CGO solutions are harmonic 
solutions defined on M" and are of the form 

Up{x) = e'''^, p e such that p ■ p = 0. 

For p = pr + ipi with pr and pi vectors in R", this means that |prP = IPiP and Pr- Pi = 0. 
When q ^ 0, CGO solutions are solutions of the following problem 

Aup + qup = 0, Up ^ e^'^ as |x| — )■ oo. (80) 

More precisely, we say that Up is a solution of the above equation with p - p = and 
the proper behavior at infinity when it is written as 

Up{x) = ef%l + ^p{x)), (81) 



27 



for ipp E Li a, weak solution of 

AV^p + 2p- V^Ap = + (82) 

The space L| for 5 G M is defined as the completion of C^(M") with respect to the norm 

II ■ II 7-2 defined as 
II iii^i- 

r - 

\\u\\l2=( {x)^^\u\^dx)\ (x) = (1 + 1x12)5. (83) 

Let — 1 < 5 < and q G L'j_^_-^ and {x)q G One of the main results in [57] is 

that there exists r] = r]{S) such that the above problem admits a unique solution with 
ipp G Lj provided that 

||(x)g||L°o + 1 < ?7|p|. 

Moreover, HV^pHia < C|p|~^||g||2,2_^^ for some C = C{6). In the analysis of many hybrid 
problems, we need smoother CGO solutions than what was recalled above. We introduce 
the spaces for s > as the completion of C^(]R") with respect to the norm || ■ \\hi 
defined as 



MHi = i / {xr\{I-Ar^u\'dx]\ (84) 

Here (/ — A)5u is defined as the inverse Fourier transform of {^)^u{^), where u{^) is the 
Fourier transform of u{x). Then we have the following 

Proposition 5.3 ([H]) Let —1<6<0 and k E W. Let q G and hence in 

—+k+e 

Hg_^_i and p be such that 

lkll^t+^+^ + l<^|p|- (85) 
Then -ipp the unique solution to (l82l) belongs to and 

|p|||^pLf+fc+-' < C\\q\\ (86) 
for a constant C that depends on 6 and rj. 

We also want to obtain estimates for ipp and Up restricted to the bounded domain X. 
We have the following result. 

Corollary 5.4 ( |18| ) Let us assume the regularity hypotheses of the previous proposi- 
tion. Then we find that 

\p\\\i^p\\H^+k+^(x) + \\^p\\h^+''+^+^{x) - ^hWn'i+'^+'ixy (8'^) 

These results show that for p sufficiently large, ipp is small compared to 1 in the class 
C'^i^X) by Sobolev imbedding. 

Let Y = H^{X) and M the ball in Y of functions with norm bounded by a fixed 
M > 0. Not only do we have that ipp is small for |p| large, but we have the following 
Lipschitz stability with respect to changes in the potential q{x): 
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Lemma 5.5 ( |16] ) Let ipp be the solution of 

and ipp he the solution of the same equation with q replaced by q, where q is defined as in 
(1271) with a replaced by a. We assume that q and q are in Ai. Then there is a constant 
C such that for all p with \p\ > \pq\, we have 

Up-^p\\Y<^\W-a\\Y. (89) 

\p\ 

This is the property used in [16] to show that a in the TAT problem ©-(IT]) solves the 
equation 

a{x) = e^f+P^-''H{x) - 'Hf[a]{x) on X, 

where 

Tif [a] (x) = a{i!f + ijf + ijfijf (x) ) , 

is a contraction map for / in an open set of illuminations; see [16]. The result in Theorem 
13.71 then follows by a Banach fixed point argument. 

5.2.2 CGO solutions and elliptic equations. 

Consider the more general elliptic equation 

- V ■ 7Vn + au = OmX, u = f on dX. (90) 
Upon defining v = ^/^u, we find that 

(A + q)v = OmX, q = - -. 

Vl 7 

In other words, we find CGO solutions for f l90|) defined on M" and of the form 



^ (a;)= 1 P-(l+^ (a;)), (91) 

a/7 



with \p\il)p{x) bounded uniformly provided that 7 and a are sufficiently smooth coeffi- 
cients. 

5.2.3 Application to qualitative properties of elliptic solutions 

Lower bound for the modulus of complex valued solutions. The above results 
show that for \p\ sufficiently large, then \up\ is uniformly bounded from below by a 
positive constant on compact domains. Note that Up is complex valued and that its real 
and imaginary parts oscillate very rapidly. Indeed, 



6' 



( cos(pj ■ x) + i sin(pj ■ x)) . 



which is rapidly increasing in the direction pr and rapidly oscillating in the direction 
Pi. Nonetheless, on a compact domain such as X, then \up\ is uniformly bounded from 
below by a positive constant. 

Let now fp = Up^g^ the trace of the CGO solution on dX. Then for / close to fp 
and u the solution to, say, f l79|) or fl90|) . we also obtain that |n| is bounded from below 
by a positive constant. Such results were used in [59] . 
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Lower bound for vector fields. For vector fields, we have the following result 

Theorem 5.6 ( |18] ) Let Up^ for j = 1,2 be CGO solutions with q as above for both pj 
and k > 1 and with Cq^\pi\ < \p2\ < co|pi| for some cq > 0. Then we have 

P ■■= :^e-^''^+'^^"Kv«p, - Up^Vup,) = ^^-^ + h, (92) 

Z\Pl\ \ / A\pi\ 

where the vector field h satisfies the constraint 

fhWc^x) < r^, (93) 
IPiI 

for some constant Cq independent of pj, j = 1, 2. 

With P2 = /Oi so that Up^ = u^, the imaginary part of (1^ is a vector field that does 
not vanish on X for |pi| sufficiently large. Moreover, let Up-^ = v + iw and Up^ = v — iw 
for V and w real-valued functions. Then the imaginary and real parts of fl92l) are given 
by 

1 O- 

3/3 = —e-^^P'-''{wVv - vVw) = ^ + Qh, 3fJ/3 = 0. 
IPiI IPiI 

Let Ui and U2 be solutions of the elliptic problem f l79|) on X such that ui + iu2 on dX 

is close to the trace of . The above result shows that 

|miVu2 — U2'Vui\ > Co > in X. 

This yields ( l24l) and the result on unique and stable reconstructions in QPAT. 

The above derivation may be generalized to the vector field Pa in (132!) with applica- 
tions in elastography. Indeed let us start from (IHT!) with p = k + i\ such that k ■ 1 = 
and k := |k| = |1|. Then using Corollary 15. 4[ we find that the following holds 



^Up = e^-^c + ip'p), Qup = e^-%5 + ip'p) 

V^Up = A;e^-^(ck-si + xp, V^Up = A;e^-^(sk + cl + Xp) 



where c = cos(l-x), s = sin(l-x), k = 1 = ^ and where |p||C| is bounded as indicated 

in Corollary [531 for C ^ {v?p, V^p, Xp, Xp}- 

Let Ml on dX be close to ^Up. Then we find by continuity that |Vmi| is close to 
|V3?Mp| so that for k sufficiently large, we find that 

|Vmi|>Co>0 inX. (95) 

This proves that H{x) = 7|Vu| is bounded from below by a positive constant provided 
that the boundary condition / is in a well-chosen open set of illuminations. 
For the application to elastography, define now 



(3a = QupV^Up — a^UpV^Up, a > 0. 

TLin{] 

4(1+q) ^""^ "5 ^ L'^P' ^P' Ap! ApJ 



Then for \p\ > pa sufficiently large so that |C| < ^'".'nl'"}^ for C e {vl, ¥^1, X'n, Xl] , we 



verify using (1M|) that 

\(3a\ > A;e2'^-"^((c5(l -«))'+ (s^ + ac^)) > ke^''-''^{mm{l,a))\ (96) 
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This provides a lower bound for /3q, uniformly on compact sets. For an open set of 
illuminations (/i, /2) close to the traces of (SjMp, ^Up) on dX, we find by continuity that 
the vector field Pa = uiWu2 — au2Vui in (132|) also has a norm bounded from below 
uniformly on X. 

Lower bound for determinants. The reconstruction in Theorem 14.71 requires that 
the determinants in (17^ be bounded from below. In specific situations, for instance 
when the conductivity is close to a given constant, such a determinant is indeed bounded 
from below by a positive constant for a large class of boundary conditions. However, it 
has been shown in [20] that the determinant of the gradients of three solutions could 
change signs on a domain with conductivities with large gradient. Unlike what happens 
in two dimensions of space, it is therefore not possible in general to show that the 
determinant of gradients of solutions has a given sign. However, using CGO solutions, 
we can be assured that on given bounded domains, the larger of two determinants is 
indeed uniformly positive for well-chosen boundary conditions. 

Let Up{x) be given by (^11) solution of the elliptic problem ( pUj) . Upon treating the 
term ipp and its derivative as in flMj) above and making them arbitrary small by choosing 
p sufficiently large, we find that ^J^Up = e^'^ + l.o.t. so that to leading order, 

^/y\/Up = e^ ^'ik + i\) ( cos(l ■ x) + i sin(l ■ x)) + l.o.t., p = k + il. 

Let n = 3 and (61,62,63) a constant orthonormal frame of M^. It remains to take the 
real and imaginary parts of the above terms and choose k = 62 or k = 63 with 1 = 61 to 
obtain, up to normalization and negligible contributions (for /c = |k| sufficiently large), 
that for 

Si = 62 cos kxi — 61 sin kxi S2 = 61 cos kxi + 62 sin kxi 

S3 = 63 cos /cxi — 61 sin /cxi 5*4 = 62 cos + 63 sin |/c|a:i, 

we verify that det (5*1, 5^2, 5*3) = — cos/cxi and that det(S^i, 5*2, 6*4) = — sin/cxi. Upon 
changing the sign of S3 or S4 if necessary to make both determinants non-negative, we 
find that the maximum of these two determinants is always bounded from below by a 
positive constant uniformly on X. This result is sufficient to prove Theorem 14. 7 j see 

m- 

Hyperbolicity of a Lorentzian metric. As a final application of CGO solutions, 
we mention the proof that a given constant vector field remains a time-like vector of a 
Lorentzian metric. This finds applications in the proof of Theorem 14.11 in [9]. 

Indeed, let k be a given direction in S""^ and p = ik + k-*- and Up = e'''^, once again 
neglecting ipp. The real and imaginary parts of Vup are such that 

g-k^.x-cjygp-x _ |k|^(a;), e-^^-'^me"-'' = \k\9^{x), (97) 
where 6{x) = kcosk-x + k-*- sink-x and 6^{x) = — ksink-x + k^ cosk-x. As usual. 

Define the Lorentzian metrics 

[)e = 26(^6 -I, i)e± =26^ - I. 
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Note that 6{x) and 6^{x) oscillate in the plane (k, k^). The given vector k thus cannot 
be a time like vector for one of the Lorentzian metrics for all x G X (unless X is a domain 
included in a thin slab). However, in the vicinity of any point xq, we can construct a 
linear combination iIj{x) = cos a6{x) +sina^-'-(x) for a G [0, 27r) such that 

k is a time-like vector for i)^ = 2ip ip — I , i.e., f),/,(k, k) = 2{iIj ■ k)^ — 1 > 0, 

uniformly for x close to xq] see [9] for more details. When 6{x) is constructed as Vu for u 
solution to (!79|) or (!90|) for boundary conditions / close to the trace of the corresponding 
CGO solution Up, then the Lorentzian metric ()^ constructed above still verifies that k 
is a time-like vector with ()^(k, k) uniformly bounded from below by a positive constant 
locally. 

6 Conclusions and perspectives 

Research in hybrid inverse problems has been very active in recent years, primarily in 
the mathematical and medical imaging communities but also in geophysical imaging, 
see e.g. [63] and references on the electro-kinetic effect. This review focused on time- 
independent equations primarily with scalar-valued solutions. We did not consider the 
body of work done in the setting of time-dependent measurements, which involves dif- 
ferent techniques than those presented here; see e.g. [S] and references. We considered 
scalar equations with the exception of the system of Maxwell's equations as it appears 
in Thermo-Acoustic Tomography. Very few results exist for systems of equations. The 
diffusion and conductivity equations considered in this review involve a scalar coefficient 
7. The reconstruction of more general tensors remains an open problem. 

Compared to boundary value inverse problems, inverse problems with internal mea- 
surements enjoy better stability estimates precisely because local information is avail- 
able. However, the derivation of such stability estimates often requires that specific, 
qualitative properties of solutions be satisfied, such as for instance the absence of criti- 
cal points. This imposes constraints on the illuminations (boundary conditions) used to 
generate the internal data that forms one of the most difficult mathematical questions 
raised by the hybrid inverse problems. 

What are the "optimal" illuminations (boundary conditions) for a given class of un- 
known parameters and how robust will the reconstructions be when such illuminations 
are modified are questions that are not fully answered. The theory of complex geo- 
metrical optics (CGO) solutions provides a useful tool to address these questions and 
construct suitable illuminations or prove their existence in several cases of interest. Nu- 
merical simulations will presumably be of great help to better understand whether such 
theoretical predictions are useful or reasonable in practice. Many numerical simulations 
performed in two dimensions of space confirm the good stability properties predicted 
by theory P [131 HSl HH [291 HH HSj • The two dimensional setting is special as we saw 
in section 15. 1[ Very few simulations have been performed in the theoretically more 
challenging case of three (or more) dimensions of space. Simulations in [JT] show very 
promising three dimensional reconstructions in the setting of diffusion coefficients that 
are close to the constant case, which is also understood theoretically since |Vm| then 
does not vanish for a large class of boundary conditions. 
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The main interest of hybrid inverse problems is that they combine high contrast with 
high resolution. This translates mathematically into good (Lipschitz or Holder) stability 
estimates. Ideally, we would like to reconstruct highly oscillatory coefficients with a 
minimal influence of the noise in the measurements. Yet, all the results presented in this 
review paper and the cited references require that the coefficients satisfy some unwanted 
smoothness properties. To focus on one example for concreteness, the reconstructions in 
Photo- Acoustic Tomography involve the solution of the transport equation fl26|) . which 
is well-posed provided that the vector field /3 is sufficiently smooth. Using theories of 
renormalization, the regularity of such vector fields can be decreased to W^'^ or to the 
BV category [ll[l9l[25]. Yet, miVm2— M2VM1 is a priori only in when 7 is arbitrary as a 
bounded coefficient [32]. The construction of CGO solutions presented in section [5l2] also 
requires sufficient smoothness of the coefficients. How such reconstructions and stability 
estimates might degrade in the presence of non-smooth coefficients is quite open. Note 
that many similar problems are also open for boundary- value inverse problems |60] . 

Finally, we have assumed in this review that the first step of the hybrid inverse 
problems had been done accurately. In practice, this may not quite always be so. PAT 
and TAT require that we solve an inverse source problem for a wave equation, which 
is a difficult problem in the presence of partial data and variable sound speed and is 
not entirely understood when realistic absorbing effects are accounted for jlOl |55] . In 
UMEIT and UMOT, we have assumed in the derivation in section |2] that standing 
plane waves could be generated. This is practically difficult to achieve and different 
(equivalent) mechanisms have been proposed [HI SI]. In transient elastography, we have 
assumed that the full (scalar) displacement could be reconstructed as a function of time 
and space. This is also sometimes an idealized approximation of what can be achieved 
in practice [H]. Finally, we have assumed knowledge of the current 7|Vm| in CDII, 
which is also difficult to acquire in practical settings as typically only the z component 
of the magnetic field can be constructed; see the recent review [SI] . The modeling of 
errors generated during the ffist step of the procedure and the influence that such errors 
may have on the reconstructions during the second step of the hybrid inverse problem 
remain active areas of research. 
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